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ABSTRACT 

We present the analysis of the Polycychc Aromatic Hydrocarbon (PAH) spectra of a sample 
of 92 typical star forming galaxies at 0.03 < z < 0.2 observed with the Spitzer IRS. We compare 
the relative strengths of PAH emission features with SDSS optical diagnostics to probe the 
relationship between PAH grain properties and star formation and AGN activity. Short-to-long 
wavelength PAH ratios, and in particular the 7.7 /im-to-11.3 /im feature ratio, are strongly 
correlated with the star formation diagnostics Z3„(4000) and Ha equivalent width, increasing 
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with younger stellar populations. This ratio also shows a significant difference between active 
and non-active galaxies, with the active galaxies exhibiting weaker 7.7 /im emission. A hard 
radiation field as measured by [OIII]/ Hfi and [NeIII]i5.Q f_im/ [Nell] 12.8 affects PAH ratios 
differently depending on whether this field results from starburst activity or an AGN. Our results 
are consistent with a picture in which larger PAH molecules grow more efficiently in richer media 
and in which smaller PAH molecules are preferentially destroyed by AGN. 

Subject headings: galaxies: active — galaxies: ISM — ISM: lines and bands — ISM: molecules — 
infrared: galaxies — techniques: spectroscopic 



1. INTRODUCTION 

The mid-infrared spectra of star forming galax- 
ies are punctuated by a series of broad peaks that 
dominate the emission between 3 and 19 fim. 
These bands are generally accepted to result 
from the vibrational modes of Polycyclic Aro- 



matic Hydrocarbons (PA Hs) ( Leger et al. 19891 : 



AUamandola et al. Il989( ) . PAH molecules — pla 



nar lattices of aromatic rings containing from lO's 
to lOO's of carbon atoms — suffuse the interstel- 
lar medium of our own Galaxy. Their vibrational 
modes are excited by absorption of UV photons, 
and in heavily star forming galaxies the PAH 
bands alone can contribute a substantial fraction 
of the reprocessed light. Given that bolometric 
infrared (IR) emission has long been an important 
measure of star formation due to the large fraction 
of dust-absorbed starlight that is re-emitted in the 
IR, the PAH bands themselves offer great promise 
as a more detailed diagnostic of star formation 
history and as a probe of star formation in heavily 
dust-obscured galaxies. 

A major factor limiting the diagnostic use of 
galactic PAH emission is that the detailed physics 
of large PAH molecules is poorly understood. 
Models of stochastic heating of dust grains pre- 
dict that the relative strengths of PAH bands 
are dependent on the size distribution of PAH 
grains and on th e ionizatio n state o f the m olecules 
(|Schutte et aTl [l993: Draine k Lil [2OO7I ). How- 
ever laboratory tests of PAHs have been limited 
to smaller size molecules than dominate the inter- 



stellar medium ( O omens et al.l .2003.: Kim et al. 



20011 ). As a result, the confidence with which the 
results of these models can be applied to interstel- 
lar PAH spectra is still uncertain. Nonetheless, 
such models indicate that measurements of rel- 
ative PAH band strengths will provide valuable 
measures of PAH growth and destruction and of 
the ambient radiation field. 



Observationally, the overall shape of the PAH 
spectrum shows remarkable similarity across a 
broad range of environments within our own 
Galaxy, and across a range of star formation 
histories in other galaxies. However there are 
still significant variations observed in the rela- 
tive strengths of certain PAH emission bands. 
Infrared Space Observatory (ISO) and Spitzer 
studies (IPeeters et all l2004t ISmith et al] [2007; 
iGalliano et al .I l2008h have revealed trends between 
PAH band ratios and various properties of the 
galaxy, including AGN activity, star formation 
history, and/or galaxy morphology. Yet, because 
these properties are themselves related, it has been 
difficult to disentangle primary physical processes 
responsible for variations in PAH spectra. 

To begin to understand the links between galax- 
ies' PAH emission and their physical properties, 
it is essential to study a sample spanning the full 
range of normal galaxy properties, and this sam- 
ple must include extensive multi-wavelength data 
to enable characterization of convolved physical 
properties. These were the primary driving goals 
behind the Spitzer SDSS GALEX Spectroscopic 
Survey (SSGSS). The SSGSS sample covers galax- 
ies from the blue cloud to the red sequence and 
transitional galaxies in between, spanning two or- 
ders of magnitude in stellar mass, colour, and dust 
attenuation. It combines both broad-wavelength- 
coverage and high-resolution IRS spectroscopy 
with a thorough suite of multi-wavelength data 
spanning the far ultraviolet to the far infrared. 
This makes it the ideal data set for studying the 
connection between galaxies' detailed IR emission 
and their physical properties. 

In this paper we study the effect of star for- 
mation, mctallicity, radiation field, and AGN in- 
cidence on PAH molecules by looking at the con- 
nection between relative PAH emission strengths 
and optical diagnostics of these properties. In Sec- 
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tion [2] we describe the sample selection, data re- 
duction and spectral decomposition. In Section [3] 
we study the relationships between the relative 
strengths of the strongest PAH bands and a range 
optical diagnostics. In Section 2] we present our 
conclusions. 

2. OBSERVATIONS AND SPECTRAL 
ANALYSIS 

2.1. Sample, Observations and Data Re- 
duction 

The full details of the SSGSS sample selection, 
observations, and data reduction are reported in 
O'Dowd et al. (2009, in preparation). To summa- 
rize: 

SSGSS is a Spitzer spectroscopic survey of 101 
normal, star-forming galaxies from the Lockman 
Hole region. This ~ 10 square degree field of low 
Galactic Hl/cirrus is extensively surveyed in mul- 
tiple wavebands. SSGSS galaxies were selected 
to have coverage by SDSS, GALEX and IRAC (4 
channel) and MIPS (24 /im, 70 /im). We applied a 
surface brightness limit of Is.g^m > 0.75 MJy sr~^, 
and a flux limit of F24p,„ > 1.5 mJy. These 
criteria yielded 154 galaxies, from which we se- 
lected 101 galaxies to span the range of normal 
galaxy properties by uniformly sampling £>„(4000) 
vs. NUV-K space. This sample has a redshift 
range of 0.03 < z < 0.2 with median redshift 
Zmed = 0.08, and a total infrared luminosity range 
of 3.7 X W^Lq < Ltir < 3.2 x 10"Lq, with me- 
dian 3.9 X 10^° Lq. 

Spectroscopy was obtained in staring mode us- 
ing the 'low-res' Short-Low (SL) and Long-Low 
(LL) IRS modules for the entire sample, provid- 
ing 5.2 — 38.0 fiui coverage with resolving power of 
~60 to 125. 'Hi-res' Short-high (SH) spectroscopy 
was obtained for a subsample of the brightest 
33 galaxies, providing 9.9 — 19.6 ^m coverage at 
R ss 600. For the bright subsample, exposure 
times were 8 min for the SL and LL modules, and 
16 min for SH. For the remaining sample, expo- 
sure times were 8 min for SL and 16 min for LL. 
All data were taken during Spitzer IRS campaign 
37. 

We used the 2-D Spitzer data products pro- 
cessed by the Spitzer Pipeline version S15.3.0, 
which performed standard IRS calibration (ramp 
fitting, dark subtraction, droop, linearity correc- 



tion and distortion correction, flat flelding, mask- 
ing and interpolation, and wavelength calibra- 
tion). Sky subtraction was performed manually, 
with sky frames constructed from the 2-D data 
frames, utilizing the shift in galaxy spectrum po- 
sition between orders to obtain clean sky regions. 
IRSCLEAN (v. 1.9) was used to clean bad and 
rogue pixels. SPICE was used to extract 1-D spec- 
tra, and these spectra were finally combined and 
stitched manually by weighted mean. 

Due to problems with some observations, the 
sample studied in this paper consists of 92 galax- 
ies observed with the low-res modules and a subset 
of 32 observed at hi-res. In addition, for subset of 
sources spectrophotometric calibration failed for 
the SL 2"'^ order module, resulting in fluxes a fac- 
tor of ~ 2 lower than expected based on the SL 
-^nd Qj-^Qj- module spectrum and IRAC photome- 
try. We eliminate these sources from analyses in- 
volving the affected spectral region, as described 
in Section [2T31 

Figure [T] shows the composite low-res spectrum 
for this sample. PAH features and complexes at 
6.2, 7.7, 8.6, 11.3, 12, 12.6, and 17 /im are promi- 
nent, as are a number of atomic emission lines. 
A thermal dust continuum component dominates 
red ward of ~ 13 fim. These composite spectra are 
normalized between 20 /im and 24 /tm. 

2.1.1. Aperture Effects 

The short- and long-low resolution IRS mod- 
ules have respective slit widths of '-~^3'.'7 and 10'.'5, 
corresponding to physical sizes of 2 — 12 kpc and 
10 — 34 kpc over the redshift range of the sam- 
ple, with respective medians 5.5 kpc and 15.6 kpc. 
The reduction process makes a flux correction for 
aperture effects assuming a point-like source. The 
maximum correction is at the intersection of the 
SL and LL modules, where the slit-loss is ~15% 
and '^45% respectively. Some error may be in- 
troduced due the fact that the different aper- 
tures sample different regions of the galaxy. For- 
tunately, the PAH features of primary interest 
are between 6.2 and 11.3 /im, and so all lie in 
the short-low regime, and are thus largely self- 
consistent. 

Aperture effects may also arise from the 
wavelength-dependent PSF width. In extended 
sources, this may result in increased sampling of 
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Fig. 1. — Composite of the 92 SSGSS low-res spectra used in this paper (black line). Blue and red lines are 
the 16th and 84th percentile of the 2.5 sigma-clipped distribution at each wavelength. Spectra are normalized 
between 20 /im and 24 fim. 



central regions at the expense of extended regions 
with increasing wavelength. We use IRAC and 
MIPS photometry in the 6, 8, 16 & 24 ^m bands 
to explore this effect, and find that the differ- 
ence due to aperture effects in the 7.7 //m-to- 
11.4 fim PAH ratio is on average less than 15%, 
and we find no systematic bias in this effect with 
any of the optical properties studied in this paper. 

The SDSS fibre diameter is 3", corresponding to 
a physical 1.6 — 9.7 kpc, with a median of 4.5 kpc, 
over the redshift range of the sample. This is suf- 
ficiently close to the scale sampled by the SL slit 
that aperture effects with SDSS data are expected 
to be minimal. 



2.2. Spectral Fitting and PAH Strengths 

To measure the strengths of the PAH features 
w e use the PAHFIT spectral decomposition code 
of lSmith et all jioO?) (S07). This code performs 
fitting of multiple spectral components, in- 
cluding PAH features modeled as Drude profiles, 
the thermal dust continuum, starlight, prominent 
emission lines, and dust extinction. For ease 
of comparison with the results of S07, we use 
identical temperatures for the thermal continuum 
components and identical central wavelengths and 
widths for the PAH features. 

We report a sample of the fitted fluxes of the 
most prominent PAH features and PAH feature 
complexes in Table [TJ These include the 6.2 and 
8.3 /.tm bands, which are discrete PAH features, 
and the 7.7, 11.3, and 17 /xm complexes, which 
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are blends of 3, 2, and 4 subfeatures respectively. 
Henceforth, we refer to individual PAH features 
and complexes of multiple features simply as 'fea- 
tures'. In Table [1] also report the fitted dust ex- 
tinction optical depth, Tg.7 and the integrated line 
strengths for [NeH]i2.8 and [NeHIJis.e ^m. The 
full table, including the entire data set and un- 
certainties for all parameters, is available in the 
Astrophysical Journal electronic edition. 

Figure [5] shows examples of PAHFIT decompo- 
sitions of the SSGSS spectra. 

2.2.1. Dust Extinction 

Dust extinction can have a pronounced effect 
on PAH feature strengths due to the strong sil- 
icate absorption features at ^10 and ^18 /im. 
PAHFIT provides the option to fit a fixed extinc- 
tion template as a free parameter. The standard 
template, which we utilize, models absorbing dust 
mixed with the emitting stars and grains, as op- 
posed to an absorbing screen. Introducing this free 
parameter has the potential to lead to degenerate 
solutions and hence to large uncertainties in PAH 
feature strengths, and so we investigated the ef- 
fect of inclusion/omission of extinction fitting on 
the fit results. We find that the strengths mea- 
sured for the 6.2 /im and 7.7 ^vci PAH features 
are very similar whether or not dust extinction is 
included, with average deviations of less than 5% 
and no systematic bias. On the other hand, the 
8.6 /im and 11.3 /xm feature strengths appear to 
be systematically overestimated without the in- 
clusion of silicate absorption by an average of 40% 
and 20% respectively. 

Due to this systematic effect, the results that 
we present in this paper utilize PAHFIT fits in- 
cluding extinction fitting. 33% of galaxies provide 
better fits with some level of silicate absorption, 
with a median fitted extinction of r = 0.7. Hav- 
ing repeated all analyses excluding silicate absorp- 
tion we find that our primary results are essen- 
tially unchanged. We make only minimal use of 
the strongly affected 8.6 /xm feature. 

2.2.2. Spline Fitting 

A commonly used alternative method for mea- 
suring PAH strengths subtracts a fitted spline to 
represent the continuum and integrate the residu- 
als within each PAH band. This method has the 



advantage of being independent of assumptions 
about the physical contributions to the fitted spec- 
trum. However, having applied this method to 
all spectra, we find that it significantly underesti- 
mates PAH strengths because it misses the PAH 
features' extended wings, and because the contin- 
uum is typically overestimated as the spline fit- 
ting points usually include substantial PAH flux. 
Worse, we find that PAH ratios are particularly 
sensitive to the systematic errors of this method 
as it depends strongly on the variable widths and 
level of blending of PAH features. As a result, we 
only report results based on the PAHFIT decom- 
positions. 

2.2.3. Problematic 6.2 fim Feature Fits 

The lower sensitivity of the 2"'' order of the 
IRS SL module results in low signal for a number 
of sources, and hence to problems in fitting the 
6.2 /im feature in some cases. Although some 
6.2 /im PAH flux is clearly present in almost all 
galaxies, PAHFIT failed to find a good fit to this 
feature in 18 cases. In many of these cases, the 
fit significantly underestimated the 6.2 /im flux, 
perhaps because of a small mismatch between the 
width or peak position of the feature and the fitted 
Drude profile. 

It is possible to estimate PAH fluxes indepen- 
dently of PAHFIT, however for the analysis of 
PAH ratios we restrict ourselves to self-consistant 
PAH flux measurements using PAHFIT. Hence, 
we omit failed 6.2 /im flts from analyses that di- 
rectly involve this feature. By BPT designation 
(see Sect. 12.2^ these include 10 non-active galax- 
ies, 3 composite sources, and 5 AGN. 

As mentioned in Section[2T] a calibration prob- 
lem with the SL 2"'' order resulted in aberrantly 
low fluxes in this order for a further subset of the 
galaxies. We identify 14 galaxies whose IRAC 
flux in the 5.7 /im band is significantly greater 
than the SL 2""^ order flux integrated across this 
band. We omit these galaxies from analyses of 
the 6.2 /xm feature. These include 10 non-active 
galaxies and 4 composite sources. 

In total, 30 galaxies are omitted from analyses 
and plots involving the 6.2 /xm feature. While this 
introduces unquantifiable biases, we still include 
such analyses where we feel that they remain in- 
formative. However the focus of this work is on the 
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Table 1 



SSGSS 


RA 


DEC 








^6.2 






^7.7 






-^8.6 




-Fll. 


i 








F 




^[Nelll] 


"^9 


.7 


1 


160 


34398 


58.89201 





066 


1 


35E- 


16 


4 


35E 


16 


8 


89E 


17 


1 


09E 


16 


4 


72E 


17 


1 


08E 


17 


1 


17E 


18 






2 


159 


86748 


58.79165 





045 


7 


04E- 


17 


2 


14E 


16 


4 


88E 


17 


5 


43E 


17 


4 


21E 


17 


3 


51E 


18 


4 


58E 


18 






11 


162 


41000 


59.58426 





047 


4 


56E- 


17 


1 


58E 


16 


3 


32E 


17 


5 


37E 


17 


3 


40E 


17 


2 


67E 


18 


2 


02E 


19 






12 


162 


36443 


59.54812 





072 


1 


55E- 


16 


6 


61E 


16 


1 


37E 


16 


1 


99E 


16 


1 


40E 


16 


3 


41E 


17 


3 


92E 


18 


1 


47 


14 


162 


52991 


59.54828 





153 


8 


67E- 


17 


3 


llE- 


16 


6 


74E 


17 


8 


93E 


17 


4 


97E 


17 


1 


06E 


17 


1 


93E 


18 






15 


161 


78737 


59.63707 





153 


1 


60E- 


17 


5 


81E- 


17 


2 


67E 


17 


3 


88E 


17 


1 


41E 


17 


8 


37E 


18 












16 


161 


48123 


59.15443 





072 


8 


38E- 


17 


3 


lOE 


16 


6 


03E 


17 


8 


49E 


17 


6 


12E 


17 


1 


35E 


17 


1 


42E 


18 






17 


161 


59111 


59,73368 





047 


5 


68E- 


16 


2 


02E 


09 


4 


HE 


16 


4 


74E 


16 


2 


81E 


16 


4 


71E 


17 


1 


HE 


17 





94 


27 


161 


11412 


59,74155 





072 


1 


32E- 


16 


1 


19E 


09 


2 


67E 


16 


3 


29E 


16 


1 


97E 


16 


6 


58E 


17 


1 


38E 


17 






28 


161 


71980 


56,25187 





103 


1 


98E- 


18 


1 


95E 


16 


4 


6SE 


17 


8 


54E 


17 


2 


67E 


17 


7 


68E 


18 


4 


23E 


18 






30 


162 


26756 


56,22390 





046 


1 


13E- 


16 


4 


02E 


16 


7 


74E 


17 


1 


OOE 


16 


3 


OlE 


17 


5 


lOE 


18 


4 


63E 


18 






32 


163 


00845 


56.55043 





117 


7 


61E- 


17 


2 


64E 


16 


5 


72E 


17 


6 


91E 


17 


1 


75 E 


17 


9 


13E 


18 


4 


76E 


18 






33 


161 


92709 


56.31395 





185 


2 


99E- 


17 


1 


61E 


16 


2 


68E 


17 


5 


48E 


17 


3 


89E 


17 


7 


02E 


18 


2 


23E 


18 





66 


34 


161 


75783 


56.30670 





046 


9 


03E- 


17 


2 


76E 


16 


6 


15E 


17 


7 


29E 


17 


3 


47E 


17 


1 


OlE 


17 


3 


14E 


18 





20 


39 


162 


04231 


56.38041 





074 


5 


47E- 


17 


1 


97E 


16 


4 


41E 


17 




54E 


17 


1 


49E 


17 


8 


42E 


18 


1 


47E 


18 






45 


161 


76901 


56.34029 





113 


2 


41E- 


17 


2 


OlE 


16 


4 


64E 


17 


6 


85 E 


17 


3 


07E 


17 


1 


73E 


17 


1 


02E 


18 






47 


163 


39658 


56.74202 





102 


7 


47E- 


17 


2 


48E 


16 


5 


14E 


17 


6 


64E 


17 


2 


89E 


17 


1 


46E 


17 


2 


18E 


18 






48 


163 


44330 


56.73859 





200 


7 


44E- 


17 


2 


92E 


16 




S6E 


17 


7 


52E 


17 


1 


88E 


17 


7 


51E 


18 


1 


98E 


18 






54 


163 


26968 


56.55812 





115 


2 


12E- 


16 


7 


50E 


16 


1 


54 E 


16 


1 


96E 


16 


1 


30E 


16 


1 


41E 


17 


1 


62E 


18 





82 


61 


163 


19810 


56.48840 





073 


3 


14E- 


17 


1 


23E 


16 


2 


OOE 


17 


4 


14E 


17 


4 


56E 


17 


3 


36E 


17 


3 


98E 


18 






62 


163 


09050 


56.50836 





133 


1 


03E- 


16 


3 


82E 


16 


7 


03 E 


17 


9 


18E 


17 


5 


97E 


17 


2 


15E 


17 












64 


163 


53931 


56.82104 





073 


2 


62E- 


16 


9 


32 E 


16 


1 


85E 


16 


2 


04E 


16 


1 


69E 


16 


4 


45E 


17 


6 


49E 


18 


1 


14 


65 


158 


22482 


58.10917 





118 


1 


82E- 


16 


6 


69E 


16 


1 


26E 


16 


1 


74E 


16 


1 


OOE 


16 


2 


72E 


17 


4 


82E 


18 





29 


69 


159 


04880 


57.72258 





076 


8 


83E- 


17 


2 


46E 


16 




92E 


17 


7 


23E 


17 


4 


54E 


17 


1 


38E 


17 


1 


91E- 


18 






70 


159 


34668 


57.52069 





090 


4 


72E- 


17 


2 


17E 


16 


4 


75E 


17 


6 


48E 


17 


2 


57E 


17 


8 


37E 


18 


2 


46 E- 


18 






73 


158 


91122 


57.59536 





080 


1 


88E- 


17 


9 


81E 


17 


2 


06E 


17 


2 


47E 


17 


2 


50E 


17 


7 


12E 


18 












77 


158 


91344 


57.71219 





044 


3 


93E- 


17 


1 


62E 


16 


3 


22E 


17 


3 


56E 


17 


1 


60E 


17 


1 


12E 


17 




46E 


18 






83 


159 


73558 


57.26361 





119 


3 


13E- 


19 


3 


06E 


16 


6 


76E 


17 


8 


36E 


17 


4 


79E 


17 


4 


69E 


18 


3 


78E 


18 






94 


159 


63510 


57.40035 





061 


7 


79E- 


17 


2 


77E 


16 


5 


33E 


17 


6 


52E 


17 


4 


08E 


17 


1 


44E 


17 




61E 


19 






95 


161 


48724 


57.45520 





115 


7 


48E- 


17 


2 


66E- 


16 


5 


84E 


17 


8 


49E 


17 


2 


52E 


17 


1 


33E 


17 


1 


04E 


19 





79 


98 


160 


29099 


56.93161 





050 


6 


40E- 


17 


7 


98E- 


16 


2 


18E 


16 


2 


64E 


16 


1 


37E 


16 


1 


81E 


17 


1 


OOE 


17 






99 


160 


30701 


57.08246 





077 


1 


62E- 


16 


5 


30E 


16 


1 


HE 


16 


1 


23E 


16 


8 


50E 


17 


2 


09E 


17 


1 


41E 


17 







Integrated fluxes in W m^^ for PAH features and complexes and for [NeII]i2.s and [jVe7//]i5.6 emission 
lines derived from PAHFIT spectral decomposition of the subsample of SSGSS galaxies with hi-res spectroscopy. 
Also reported are the redshift of the galaxy, z, and dust extinction given as rg.y, the optical depth at 9.7 fiui. PAH 
strengths and absorption come from fitting of low-res data while emission lines come from fitting of hi-res data. 
Results for the full data set, including uncertainties for all values, is included in the electronic edition of the 

Astrophysical Journal. 
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10 15 20 25 30 10 15 20 25 30 

Wovelength (/^m) 



Fig. 2. — Examples of PAHFIT decompositions of SSGSS spectra (black squares) in relative pF^, showing 
a range of spectral shapes. The best fit (green) is composed of thermal dust continua at set temperatures 
(red), PAH features (blue), stellar light (magenta), emission lines (purple), and silicate absorption (dotted). 
SSGSS catalog number is also given. 
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more robust PAH measurements - in particular th 
7.7 /im and 11.3 /im features. 

2.2-4-. The Relative Dominance of Star- Formation 
and AGN 

Throughout this paper we colour-code to indi- 
cate location in log [01 1 1] 5007 /H[3 v s. log \N I Ha 



space: the so-called BPT diagram of lBaldwin et al. 



(|l98l[ ) (Fig. [31). Red p oints indicate more p ower- 
ful AGN satisfying the iKewlev et all (|200l[) des- 
ignation, green points indicate composite sources 
with weaker active components betw een the Kew- 
ley and the iKauffmann et al.l (|2003l ) designation, 
while black points indicate galaxies dominated by 
stellar light, below the Kauffmann designation. 

BPT designation must be interpreted care- 
fully. Rather than an absolute discriminator be- 
tween star-forming galaxies, composite sources, 
and AGN, for this study we instead interpret it 
as a measure of the relative dominance of star- 
formation and AGN with the Sloan aperture. The 
similarity of the Sloan and Spitzer SL apertures 
means that BPT designation measures this rela- 
tive dominance within roughly the same galactic 
region responsible for the observed PAH emission. 
The regions probed in this study are larger than 
those probed in earlier, local studies. This is an 
important difference as it means that we observe 
a more global stellar population, where any AGN 
component is less overwhelming. 

2.2.5. Correlation of PAH Ratios and Compari- 
son to Other Studies 

Ratios of the integrated fluxes of the fitted 
Drude profiles were used to determine PAH lu- 
minosity ratios, designated L\i/L\2, where A is 
the central wavelength of the feature in microm- 
eters. Figure [3] shows the tight correlations be- 
tween different short-to-long wavelength PAH ra- 
tios, ie. 2/^11.3 versus L7.7/L11.3 follows a tight 
locus with ~0.25 dex scatter. Lg.e/iii.a versus 
L7.7/L11.3 shows a similar scatter, but broadens 
towards low ratio values. The tight correlations 
suggest that lower-wavelength PAH features are 
strongly coupled, and are less-strongly tied to the 
11.3 /xm feature. 

In this figure, we also compare our results to 
two other spectroscopic studies of PAH intensities: 

S07 present spectral decomposition of IRS ob- 




[NII]/Ha 



Fig. 3.— BPT plot ([Baldwin et al. Ill98lh for the 



SSGSS sample, coloured coded according to BPT 
designation; blue points are galaxies dominated 
by star formation, green points are galaxies with 
significant stellar and AGN components, and red 
points are AGN-dominated. Strictly, these desig- 
nations are measures of the relative dominance of 
AGN and star-forming components within the ob- 
served aperture. However, for convenience, we uti- 
lize the standard naming convention and refer to 
these as star-forming galaxies, composite sources, 
and AGN, respectively. The upper an d lower dot- 
ted lines represent the [Kewlev et al. 1 (l200l[) and 
Kauffmann et al.[ ( 20031 ) divides respectively. 



servations of 59 SINGS galaxies, and we use an 
identical spectral decomposition technique (see 
Sect. 12. 2p . An important difference between S07 
and our survey is, with their local targets and 
shorter exposure times, they are dominated by 
the highest surface brightness regions near the 
nucleus, while our observations sample more of 
the extended structure. Their reported PAH ra- 
tios are very similar to ours, with the exception 
that the S07 sample have, on average, higher 
L7.7/L11.3, and lack the small number of galax- 
ies with Le. 2/^11. 3 < 0.3 that we observe. 

Galliano et al] (|2008l ) (G08) present spectral 



decomposition of ISO and Spitzer IRS observa- 
tions of Galactic HII regions and dwarf spiral 
and starburst galaxies. We take the results of 
their most comparable decomposition technique. 



8 



0.2 



1.0 

^.s/^l 1.3 



0.2 1.0 

'^B.e/^ii.a 



Fig. 4. — The PAH band ratios Le. 2/^11 3 versus L7 7/L11 3 (left) and Lg.e/iii.a versus L7. 7/^11.3 (right). 
Yellow triangles are SINGS galaxies from Smith et al.l ( 2007 ). derived using PAHFIT with the same param- 
eters. The dashed and dotted lines show the best linear regression and the dispersion for the PAH ratios 
of Galactic HH regions and dwarf spiral and starburst galaxies found bv iGalliano et al" ( 2008f ). fit using a 
different technique. Colour-coding is the same as in Figure [31 



in which fit they Lorenzian profiles to the PAH 
bands (as opposed to the Drude profiles used by 
PAHFIT). While their sample show very simi- 
lar trends in the 7.7 /im and 8.6 /im features 
versus the 11.3 fim feature, they find that the 
6.2 /im feature is on average stronger by a factor 
of 1.3 to 2 than both our results and those of S07. 
This discrepancy is likely due to the differences in 
the fitting techniques. 

3. RESULTS AND DISCUSSION 

3.1. PAH Ratios as Measures of Grain 
Size and Ionization 

Models of the stochastic heating of dust grains 
(Tielend 120051: iDraine &: Lil l2007t Ischutte et al. ' 



1993|) show that the relative power emitted in 
a given PAH band is strongly dependent on the 
distribution of grain sizes. In general, smaller 
dust grains emit more power in the shorter wave- 
length bands and larger grains dominate longer 
wavelengths bands. In the regime of grain sizes 
expected in interstellar dust (lO's to lOO's of 



carbon atoms), rapid drops are expected in the 
6.2 /im and 7.7 /im bands relative to longer- 
wavelength features, and in 6.2 /im relative to 
7.7 /tm, with increasing grain size. 

It is also expected that the ionization state 
of a PAH molecule will have a dramatic effect 



on its spectrum (|Draine fc Lil 120071 ). In particu- 
lar, carbon-carbon vibrational modes are known 
to be signi ficantly more intense in ionized PAH 
molecules (|Tielens|[2005l) . As the 6.2 /im and 
7.7 /im bands result from radiative relaxation of 
CC stretching modes, the ratios of these bands to 
those arising from carbon-hydrogen modes such 
as the 11.3 /tm feature are expected to drop by 
an order of magnitude between completely neu- 
tral and completely ionized PAH clouds. At the 
same time, the 6.2 /im feature should change lit- 
tle relative to the 7.7 fim feature as the ionization 
fraction changes. 

Comparison of the PAH band power ratios 
^6. 2/^7. 7 to L11.3/L7.7 is useful in extricating 
these two effects, with the former being sensitive 
to PAH grain size but relatively unaffected by ion- 



ization as both bands result from the similar CC 
modes, while the latter is very sensitive to ioniza- 
tion, and somewhat less so to grain size. Figure [5] 
shows these r atios for our sample plotted over the 
Draine fc Lil (j2007r ) models. Most of the sample 
is restricted to a tight locus between 0.2 and 0.4 
in both ^6.2/^7.7 and in. 3/^7. 7, although within 
this locus there is a weak trend in the direction ex- 
pected for constant ionization fraction with chang- 
ing grain size, albeit with significant scatter. We 
discuss the distribution of active galaxies in this 
plot in Section [33] 



3.2. 



Comparison of PAH Ratios with Op- 
tical Diagnostics 



We compare the PAH spectra resulting from 
our PAHFIT models to SDSS optical diagnos- 
tics of star formation history and AGN activ- 
ity. Figure [5] shows the ratio of short-to-long 
wavelength PAH bands (^5.2/^7.7, ^e. 2/^11. 3, 
and Lt,i/Lii,3) versus the star formation diag- 
nostics I?ri(4000), Specific Star Formation Rate 
(SSFR), a nd Ha Equivalent Wid th (Ha EW). 



i:)„(4000) (jKauffmann et al.l l2003l ) measures the 
4OOOA break over a narrower bandpas s than the 



standard 1^(4000) (|Bruzual et al.lll983l ). and pro- 
vides a meas ure of the luminosity- wei ghted stellar 
age. SSFR (jBrinchmann et all \2004) is the star 
formation rate as determined by population syn- 
thesis models per unit stellar mass. The Ha line 
provides a direct measure of recent star formation, 
and in particular the presence of short-lived OB 
stars. 

These star formation diagnostics are relatively 
free of the effects of dust extinction; Z?„(4000) 
and Ha EW because of the narrow bandpasses 
over which these diagnostics are calculated, and 
SSFR because it has been corrected for extinction 



(jBrinchmann et al.ll2004[ ) 



The most striking correlations are seen rela- 
tive to L7. 7/^11.3, with the short wavelength PAH 
band becoming more dominant with increasing re- 
cent star formation. The ratios of the 6.2 fim fea- 
ture to longer wavelength features follow similar, 
albeit weaker trends to L7.7/L11.3. 

Table [2] shows the results of the regression anal- 
yses of PAH ratios to optical properties. We use 
Kendall's correlation to determine the significance 
of the trends, given as probabilities of obtaining 



the observed data given the null hypothesis of no 
correlation. We obtain a significance of <5% for 
L6.2/Lii,3 and L-j,t/ Ln,^ versus all star formation 
diagnostics, while ^6.2/^7.7 shows no significant 
correlations. 

[0III]/H(3 provides a measure of the ioniza- 
tion parameter weighted by SFR, and hence of 
the hardness of the radiation field present in the 
galaxy (Kewlev et al..,2001.) . From Figure [7] (left) 
it can be seen that there is a global trend between 
this line ratio and ^7. 7/^11. 3, and no strong trends 
with ratios involving the 6.2 /im feature. 

The emission line ratio [Nil]/ Ha provides a 
diagnostic of gas-phase metallicity, although it is 
also affected by radiation field hardness. This ra- 
tio strongly correlates with the age of the stellar 
population. As expected, [Nil]/ Ha follows the 
trends observed with star formation diagnostics 
(Fig. [71 middle). The best global correlation is 
again with L7.7/L11.3, with the ratios of short-to- 
long wavelength PAH bands decreasing with in- 
creasing emission line ratio. 

Table [2] also shows the results of the correlation 
analysis for PAH ratios versus emission line ratios. 

A more reliable measure of gas-phase metal- 
licity, der ived from a range S DSS optical nebu- 



lar lines ( Tremonti et al.l [20041 ) is available for a 



subset of the sample, although contain no BPT- 
designated AGN or composite sources (Fig. [71 
right). The correlations of PAH ratios with this 
metallicity measure are slightly weaker than that 
for [Nil]/ Ha, suggesting that the trends ob- 
served in [NII]/Ha may be dominated by the 
trend with radiation hardness. However there is 
a notable locus of high short-to-long wavelength 
PAH ratios for lower mctallicitics. Galaxies in 
these upper left loci have the lowest metallicities 
(12 + log{0/H) < 9) and are among the youngest 
in the sample (Z?„(4000) < 1.3), and exhibit the 
hardest radiation fields of the star-forming galax- 
ies {[Oil I] /HP > 0.4). 

It is known that PAH emission aro und 8 /im is 

stron ger in higher metallicity galaxies (jEngelbracht et al, 
20051 ). The dominance of short wavelength PAH 



features in the lowest metallicity galaxies as seen 
in Figure [7] (right) may result from a difficulty 
in building large PAH molecules in an under- 
rich medium. It may also result from other 
age-dependent effects given the close link be- 
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Fig. 5. — The PAH band ratios £6.2/-^7.7 versus L11.3/L7.7. Grey lines show the expected PAH ratios as a 
function of grain size and ionization state from the models of iDraine fc Lil (j2007l ). The lines represent the 
expected ratios for fully neutral or fully ionized PAH molecules of a given number of carbon atoms. These 
are meant to be illustrative only, as real PAH ensembles will contain mixtures of ionization states and grain 
sizes. Colour-coding is the same as in Figure [3] 



tween metallicity and age. However it is not 
expected to arise through the link between metal- 
licity and radiation field hardness; a hard radia- 
tion field will either preferentially destroy small 
PAH molecules, decreasing this ratio, or ionize 
molecules, which should not produce the observed 
change in i6.2/-^7.7 with metallicity, as seen in 
Section 13.11 

3.3. Separating the EfTects of AGN Activ- 
ity and Star Formation 

From Figures [6] and [7] it can be seen that galax- 
ies with dominant AGN components (those above 
the Kauffmann divide on the BPT diagram) have. 



on average, weaker short-to-long wavelength PAH 
bands. There is a strong link between AGN inci- 
dence and the star formation history of a galaxy; 
most notably, AGN are not apparent in the most 
strongly star- forming galaxies, because of swamp- 
ing of diagnostic lines, dust obscuration, or for 
reasons related to the life-cycle of the AGN. The 
link between age and BPT designation is apparent 
from Figure[6l This makes it difficult to determine 
whether AGN activity, star formation history, or 
a combination of the two are responsible for the 
trends observed between PAH ratios and optical 
properties. 
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Fig. 6. — PAH ratio versus -D„(4000), SSFR, and Ha equivalent width. Colour-coding is the same as in Fig- 
ure [3l The lines show the median values for galaxies with no AGN component (dashed) and AGN/composite 
sources (dotted). 



3.3.1. The Effect of Star Formation 

We repeat the regression analysis and linear 
fits from Section 13.21 independently for subsam- 
ples of star-forming galaxies and AGN/composite 
sources. For both subsamples, the correlations 
with L7.7/L11.3 still hold for Ha EW. However, 
for Z?„(4000) and SSFR the correlation remains 
significant at better than the 5% level only for the 
star-forming subsample. For Lq ^/Lh,^ the corre- 
lation only holds for the star formation-dominated 
subsample with respect to Ha EW. 

In general, it appears that star formation his- 
tory does have an influence on ratio of the 7.7/im- 
to-11.3 ixm PAH features, independent of an AGN 



component. Table [2] includes the regression anal- 
ysis for the split samples. 

3.3.2. The Effect of AGN 

The split-sample regression analysis for Ly y/Ln a 
vs. [OIII]/ H(3 reveals that the correlation still 
holds for the AGN-|-composite subsample with a 
similar slope to the full sample. While there is a 
formal correlation for the star-forming subsample, 
the best-fit slope is so shallow as to suggest that 
the correlation is not meaningful. This suggests 
that there is a correlation between the hardness of 
the radiation field originating from an AGN and 
the 7.7 /Ltm-to-11.3 /zm ratio. 

In Figure [8] (left) we show our sample on the 
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Fig. 7. — PAH ratio versus [Nil]/ Ha and [OIII]/ H p. Colour-coding is the same as in Figure [3l The 
Hues show the median values for galaxies with no AGN component (dashed) and AGN/composite sources 
(dotted). 



BPT diagram, colour-coded as in Figure [31 For 
better visual clarity, symbol size is proportional 
to the ratio of long-to-short wavelength PAH band 
luminosity (upper: Ly. 7/^6.2, lower: in. 3/^7. 7), 
and hence to average grain size and/or neutral 
PAH fraction. The histograms (Fig. [51 right) show 
the distribution of PAH ratios for the different 
BPT designations. The strongest trend is found 
for L7.7/L11.3, for which the active galaxies show 
significantly lower values than star-forming galax- 
ies (with <0.1% significance by a KS test). More 
dominant AGN satisfying the Kewley designation 
have only marginally lower L7.7/L11.3 than com- 
posite sources satisfying the Kauffmann designa- 
tion (10% significance). For ^6.27-^7. 7, the differ- 



ences in the subsamples are not highly significant; 
there is a 10% significance for the difference in this 
ratio between active and star-forming galaxies. 

Although there is only a tentative difference 
in the L-i/j/Ln^^ distributions between the BPT- 
designated full AGN and composite sources, these 
galaxies have indistinguishable distributions of 
i:>„(4000) and Ha EW by a KS test. Thus, any 
difference in PAH ratios is likely due to AGN dom- 
inance and not to a link between AGN and star 
formation properties. 

To better establish the independence of this 
relationship to star formation we analyse a con- 
trolled subsample restricted to a narrow age range. 
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Diagnostic Sample 



ifi 



A. (4000) 



Ha EW 



[Nil] 
Ha 



[QUI] 

H(3 



full sample 
star-forming 
AGN /composite 



full sample 
SSFR star-forming 

AGN /composite 



full sample 
star-forming 
AGN /composite 



full sample 

star-forming 
AGN /composite 



full sample 
star-forming 
AGN /composite 



numbers given are the slopes of the best linear fits 



1.3 



^7.7 
1.3 



-0.07 (1.6%) -0.26 (<0.01%) 
-0.39 (3%) 



0.2 (0.17%) 0.21 (<0.01%) 
0.13 (0.2%) 



0.16 (0.01%) 0.19 (<0.01%) 
0.1 (0.5%) 0.16 (0.01%) 
0.19 (2%) 



-0.5 (0.1%) -0.53 (<0.01%) 



-0.20 (1.2%) 

-0.07 (4%) 
-0.25 (5%) 



Table 2: Results of linear fits and regression analysis for PAH ratios versus SDSS optical properties. The 



followed in brackets by the significance of the trend 



according to Kendall's rank correlation. Dashes indicate that the correlation is not significant at the 5% 
level or better. 



with 1.3 < A>(4000) < 1.5. Figure [i shows the 
histogram of ^7. 7/^/^11. 3 for galaxies in this range, 
split between galaxies with dominant AGN activ- 
ity (either composite or full AGN by BPT desig- 
nation), versus galaxies with dominated by star 
formation. In this range, both subsamples have 
indistinguishable distributions of both Z3„(4000) 
and Ha EW. A KS test of these PAH ratios shows 
that they are different, with a significance level of 
2%. It is therefore likely that AGN play some role 
in reducing the relative 7.7 /^m-to-11.3 /im PAH 
emission, and that the trend is not solely due to 
the link between AGN status and the age of the 
stellar population. 

3.4. The Effect of Age, Metallicity and 

Hard Radiation Fields on PAH Molecules 

In order to better separate the role of AGN 
and star formation, we look at the relationship 
between hardness of the radiation field as mea- 
sured by [OIII]/H/3 with stellar population age 
as measured by _D„(4000) (Fig. [T0|). As a sensi- 
tive age diagnostic, I?„(4000) is highly correlated 
with [Nil]) /Ha, and so this plot is analogous to 
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Fig. 9. — Histogram of L7.7/L11. 3 for star-forming 
galaxies (blue) and sources with any AGN com- 
ponent (red), only including sources with 1.3 < 
£'„(4000) < 1.5. A KS test indicates that samples 
are different at a 2% significance level. 

the BPT plot, but with a clearer separation be- 
tween young and older stellar populations. In this 
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Fig. 8. — Left: BPT plot for SSGSS sample. Colour-coding is the same as in Figure[3l Symbol size increases 
with decreasing short-to-long wavelength PAH ra tios; iv6.2/-^7.7 (upper) a nd L7.7/L11.3 (lower). Th e upper 
and lower dotted lines in each plot represent the iKewlev et al.! ((2OOI. ) and iKauffmann et ah I (|2003l) divides 



respectively. Right: Histogram of PAH ratios (upper: ^6.2/^7.7, lower: L-j t/Lh,^) for star-forming galaxies, 
composite sources, and AGN. 



plot symbol size is again proportional to long-to- 
short wavelength PAH ratio for the purpose of 
visual clarity. Here we see that, for hard radiation 
fields as measured by [OH I] / HP, there is a strik- 
ing difference in PAH ratio between younger and 
older stellar populations. Younger stellar pop- 
ulations (D„(4000) < 1.3) with hard radiation 
fields ([OIIIJ/H/3 > 0.3) have stronger short wave- 
length PAH bands than older galaxies with simi- 
lar [0III]/H(3 (significance: 0.05% for L7.7/L11.3, 
3% for Lq,2/L7,j); they also have stronger short 
wavelength bands than galaxies of all ages with 
softer radiation fields ([OIIIJ/H/3 < 0.3, signif- 



icance: 2% for L7.7/L11.3, 1% for ig. 2/^7.7)- 
Older populations with hard ionizing fields (which 
are also predominantly AGN-dominated) have 
stronger long wavelength bands than the other 
populations (significance: <0.01% for Lt^j/Lh,^, 
8% for ^6.27-^7. 7)- The differences are illustrated 
in the histograms to the right of Figure [TOl 

Thus, it seems that the behaviour of PAH 
molecules in the presence of a hard radiation field 
differs depending on the age of the stellar popula- 
tion, and so it may depend on whether the radia- 
tion comes predominately from starburst activity 
or an AGN. These trends may also arise from a 
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Fig. 10. — Left: [OIII]/ H (3 vs. D„(4000). Colour-coding is the same as in Figure [3] Symbol size increases 
with decreasing short-to-long wavelength PAH ratios; Lq ^/L'!,^ (upper) and Lj^r/Lu^s (lower). Dashed 
lines divide sample into the subpopulations used in the histogram. Right: Histogram of PAH ratios (upper: 
Le.2/LT,7, lower: L7.7/L11.3) for both old (Z3„(4000) > 1.3) and young (Z?„(4000) < 1.3) sources with hard 
radiation fields {[OIII]/ H/3 > 0.3) and all sources with softer radiation fields {[OIII]/ H/3 < 0.3). 



link between metallicity and both radiation field 
hardness and stellar population. Below, we dis- 
cuss the factors that drive the evolution of short- 
to-long wavelength PAH ratios: 

The youngest galaxies begin their evolution 
with low D„(4000), low metallicity, a moderately 
hard radiation field resulting from starburst activ- 
ity, and high short-to-long wavelength PAH ratios. 
The increase in metallicity with age is expected 
toallow the larger grains that dominate the longer 
wavelength feature in each ratio to be assembled 
more efhciently, resulting in a drop in this ratio as 
seen in Section [ 



Additionally, as age increases the starburst ra- 
diation field also softens. This may result in a 
decrease in the ionization fraction of the PAH 
molecules, which, according to the iDraine fc Li 



([20071) models (Fig. [5|) , may contribute to the drop 
in 7.7 /zm-to-11.3 fim ratio between the youngest 
galaxies in this sample and older galaxies. 

Even after these youngest, starburst stages, ag- 
ing galaxies continue to exhibit a decrease in short- 
to-long wavelength PAH ratios, and in particular 
in L7.7/L11.3 (see Fig. |6]). Some of these galaxies 
experience AGN activity, and so are subject to ra- 
diation fields that are often much harder than in 
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Fig. 11. — [A''e7//]i5.6 tj.m/[NeII]i2.s /xm vs. L7.7/L11.3. Circle s are SSGSS galaxi es, with colour-coding is 



the same as in Figure |31 Diamonds are SINGS galaxies from ISmith et all ( 20071 ). with the red and blue 
points representing AGN/composite sources and star- forming galaxies respectively. 



star burst galaxies. If the initial decrease in PAH 
ratio with age is due to increasing neutral frac- 
tion of PAH molecules, then AGN do not seem 
to reionize these molecules; indeed, the trend of 
diminishing PAH ratio with galaxy age continues 
through an increasing incidence of AGN. 

In Section [3.3.2l we saw that AGN seem to cause 
a decrease in PAH ratio, particularly iy.y/Ln.a, 
independently of star formation diagnostics, and 
that -L7.7/L11.3 also drops with AGN-sourced ra- 
diation hardness (Fig. [71 lower left). This may 
be due to preferential destruction of small PAH 
grains by X-rays and/or shocks from the AGN. It 
may also be due to heating of PAH molecules by 



the AGN, as ionization fra ction decrease s with in- 
creasing gas temperature (lTielendl2005l) . Thus, if 
the AGN is capable of heating PAH to high tem- 
peratures in dense regions, compared to the more 
diffuse heating by stars, then the net effect may 
be to reduce the ionization fraction. 

3.4.1. [Nell] and [Nelll] Measurements 

This differing effect of an AGN- or starburst- 
sourced hard radiation field on the PAH spec- 
trum is consistent with the result of S07, who plot 
Lt.t/Lii,^ versus [A^e///]i5.6 ,™/[A^e//]i2.8 
for their sample. [NeIII]/[NeII] provides a NIR 
measure of the hardness of the radiation field. 
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For a subsample of galaxies with significant AGN 
components (Seyferts and LINERs), they find a 
strong anti-correlation between the two ratios: the 
harder the ionizing field in the AGN, the weaker 
the short wavelength PAH band. However the J 
hardness of the radiation field in star-forming j' 
galaxies — radiation that is primarily stellar in 
origin — shows no correlation with PAH ratio. 
Brandl et al. I (|20n6l ) find a similar absense of cor- 
relation for 22 starburst nuclei over a similarly 
broad range of [NeIII]/[NeII]. 

We have measured [Nell I] j [iVe//] for the sub- 
sample of 32 galaxies with the IRS hi-res spec- 
trograph. In Figure [11] we compare our results 
with S07. Our results are consistent, however the 
SSGSS sample does not include [iVe///]/[iVe//] 
values as high as those found by S07, and so 
we don't see the full trend between these ratios. 
The likely reason for the difference in measured 
[iVe///] / [-/Ve//] is that for the local galaxies stud- 
ied by S07, the IRS slit samples a smaller, more 
central, and hence more AGN-dominated region 
than for our more distant galaxies. 
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3.5. The 17 /im Feature 

The 17 /J,m PAH band is well separated 
from the denser blueward bands, and hence 
provides one of the cleanest measures of PAH 
stre ngth. This band m ay arise from C-C-C bend- 
ing (jPraine fc Lil 120071 ) . however the true vibra- 
tional mode or modes responsible for this emission 
are not established. Nonetheless, it is expected 
that the dominant contributing grain size will be 
larger than for lower wavelength bands. This band 
falls in the long-low IRS module, which has a sig- 
nificantly larger slit width than both the short-low 
module and the SDSS aperture, at lO'.'S vs. ~3'.'7 
and 3" respectively. As a result, the sampling of 
different regions of the galaxy introduces unquan- 
tifiable uncertainties, and so we limit the analysis 
of this feature to qualitative comparisons. 

In FigurefT2lwe plot the ratio of the 7.7 /im bands 
to the 17 /im band versus SDSS optical properties. 
With respect to most optical diagnostics, L7.7/L17 
displays similar correlations to L7.7/L11.3, with 
similar dispersions. This ratio shows a signif- 
icantly better correlation with respect to SSFR, 
and a worse correlation with respect to [iV//]/i/a. 
As with L7.7/L17, no correlation is observed with 
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Fig. 12. — Lu/Ln vs. SDSS optical diagnostics 
and Nelll/Nell. Colour-coding is the same as 
in Figure O 

respect to NeIII]/[NeII]. Dividing the sample 
between BPT-designated AGN/composite sources 
and star-forming galaxies, the same correlations 
are observed with all star formation diagnostics 
for star- forming galaxies. The correlations are not 
significant for other diagnostics or for the subsam- 
ple of galaxies with AGN components. 

In Figure [13] we compare the histograms of the 
PAH ratios ie. 2/^7.7, -^7.7/^11.3, -^7.7/^17, and 
Lii.s/Ln divided into subsamples of varying of 
AGN activity. As was seen in Figure [8] ^6.2/^7.7 
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Fig. 13. — PAH ratio by BPT designation, with 
grey for star- forming galaxies, green for composite 
sources and red for AGN. Top left: ^6.2/^7.7, top 
right: i7.7/-^ii.3, bottom left: ig. 2/^17- bottom 
right: L7.7/L17. 

reveals a slight relative increase in the strength 
of the longer-wavelength PAH feature for galaxies 
with an AGN component while Ly^t/Lh^^ shows a 
much stronger increase. There is also a significant 
increase in the strength of the 17 /xm feature rel- 
ative to the 7.7 ^m feature with increasing AGN 
incidence and power. There is an increase in the 
17 fim feature relative to the 11.3 fim feature 
between star-forming galaxies and full AGN, but 
not between star-forming galaxies and composite 
sources. 

This may indicate a threshold in AGN power 
needed to destroy the larger PAH grains that dom- 
inate the 11.3 /zm feature; only full AGN are ca- 
pable of destroying such grains with enough ef- 
ficiency to register a difference in the L11.3/L17 
ratio. Conversely, the average molecule sizes that 
dominate both the 6.2 and 7.7 /zm features are 
destroyed with similar efficiency even for compos- 
ite sources, resulting in little change in ^6.2/^7.7 
with AGN status. 

In Figure [TJ] we study the 17 /zm feature 



in comparison to shorter wavelength features. 
^6. 2/^7. 7 vs. L7.7/L17 shows a very narrow range 
of the former ratio over an order of magnitude vari- 
ation in the latter. Ljj/Lu correlates well with 
^7.7/^11.3- Interestingly, an even tighter correla- 
tion is observed between L-j^-j/Ln and L11.3/L17, 
suggesting that the 7.7 /im and 11.3 /im fea- 
tures are more closely tied than the 11.3 /im and 
17 jjum features. 

4. CONCLUSIONS 

We have studied the relationships between 
a range of optical diagnostics and the relative 
strengths of PAH emission bands for the SS- 
GSS sample of 92 GALEX-selecied normal, star- 
forming galaxies observed with the IRS lo-res and 
hi-res spectrographs. The observed PAH spec- 
tra a£e_consi^ent_with SINGS galaxies studied 
by|SmithetaD(|2QQ3) and largely consistent with 
the galactic H I I regi ons and galaxies studied by 
iGalliano et al. I (|2008h . although in the case of the 
latter, the different analysis methods appears to 
yield a ^^20-50% difference in the strength of the 
6.2 /im feature. 

Short-to-long wavelength PAH feature ratios 
exhibit a number of trends with the SDSS optical 
diagnostics of the galaxies' star formation histo- 
ries, metallicities, and radiation fields. The most 
striking of these is the correlation of the 7.7 /im- 
to-11.3 /im feature ratio with the star formation 
diagnostics Ha EW and D„(4000), and with the 
emission line ratios [Nil]/ Ha and [OIII]/H/3, 
and with metallicity. The correlation of this PAH 
ratio with star formation diagnostics is indepen- 
dent of the presence of an AGN component, in- 
dicating that the stellar population plays an im- 
portant role in determining the relative feature 
strengths. This is consistent with an increase in 
production of large-grain PAH molecules with in- 
creasing metallicity, and with a correlation be- 
tween the fraction of PAH molecules that are ion- 
ized and the hardness of any starburst radiation 
field. 

The presence of an AGN component, as deter- 
mined by the galaxy's location on the BPT dia- 
gram, is correlated with a reduction in the ratio 
of the 7.7 /im to 11.3 and 17 fim features. For a 
subsample with matched £'„(4000) and Ha EW, 
galaxies with any AGN component (AGN or com- 
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Fig. 14. — ^6. 2/^7. 7 (left), -L7.7/L11.3 (middle), and Lii.3/ii7 (right) vs. L7.7/L17. Colour-coding is the 
same as in Figure [31 



posite sources) have weaker relative 7.7 /im emis- 
sion than quiescent galaxies. This is consistent 
with a picture in which smaller PAH grains are 
preferentially destroyed by shocks and/or x-rays 
from the AGN. 
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